The globule-coil transition in a mean field approach 
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In this paper we study the unfolding transition observed in polymer stretching 
experiments. We use a known model and extend it with a mean-field-like interaction. 
Expressions for all necessary thermodynamic variables as a function of the model 
parameters are derived in closed form. To obtain the exact force-extension relation 
at constant temperature, we invert these relations numerically. Below a critical 
temperature the model exhibits a sharp first order phase transition. 

PACS numbers: 64.60. Cn 82.35.Lr 



I. INTRODUCTION 

The unfolding transition of polymers is observed in different types of experiments. The 
unfolding is a single- or multi-step process, depending on the experimental conditions. 

In force-clamp experiments, one is able to study the mechanical unfolding of polymers 
at relatively constant force. In such experiments, one applies a sudden force (the time to 
complete the force step is of the order of 10ms 111). With feedback techniques the force is 

n n 

then kept constant. Both experiments [1| and numerical simulations 2] show that a certain 
protein (called ubiquitin) unfolds in a single step, while an other protein (called integrin) 
unfolds in multiple steps In J], Q], force-extension relations of single DNA molecules are 
measured in the fixed stretch ensemble. In these experiments, one keeps the extension of the 
DNA molecule relatively constant, and the average applied force is measured. Depending 
on the solvent conditions force plateaus or stick-release patterns are observed. 

The unfolding of polymers is also studied under more dynamical conditions. The two 
most common are the mechanical stretching at a force that increases linearly with time 
(force-ramp) |5j and the deformation of polymers in a flow . Numerical simulations 
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FIG. 1: (color online). Plot of the force-extension relation at constant temperature (h(3 = 0.7). 
The black solid line is the stable solution. The blue dotted line corresponds with the unstable 
solution. The length of the chain is 10. 



show that under these conditions, the unfolding is mostly a multi-step process. 

The unfolding transition from compact globule to stretched coil was already theoretically 
predicted 15 years ago [9], based on heuristic arguments. Nowadays, the theoretical study 
of this transition is dominated by numerical simulations. Instead of performing numerical 
simulations, we study a simplified model that can be solved in closed form almost completely. 
Despite the simplicity of the model, it still exhibits a first order phase transition for every 
chain length at low enough temperatures. For this reason, the model is also interesting on 
a purely theoretical ground. 



In 



10 



111 ] we introduce a model to describe polymers. In the present paper, we start 



with shortly reviewing the results obtained in these papers. We proceed with deriving 
expressions in closed form for the temperature and the external force by minimizing the free 
energy. These results are used to study the unfolding transition of the polymer. The last 
part of the paper gives a discussion of the results. 



II. MODEL 



Our model is a one-dimensional persistent random walk with an external force. It is 
described in detail in References lOj and ll|. The model depends on two parameters e 
and fi, the probabilities to go straight on when walking to the right, respectively to the left. 
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FIG. 2: (color online). Plot of the average number of kinks as a function of the extension at 
constant temperature (h/3 = 0.7). The black solid line is the stable solution. The blue dotted line 
corresponds with the unstable solution. The length of the chain is 10. 



This is not a Markov chain since the walk remembers the direction it comes from. Two 
important variables are the end-to-end distance x and the number of reversals of direction k 



(number of kinks). In 



10] , the generating function method is used to obtain exact analytical 



expressions for the average end-to-end distance and the average number of kinks 
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with a the lattice parameter and n the total number of steps. Also the average of the squared 
end-to-end distance can be calculated with the same method (for the result, see appendix 
lAl) . Starting from the definition of Boltzmann, the following expression for the entropy S is 



obtained in 



ill 



S = -<*>]»£. 
2 
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(We use units for which k B = 1.) 



III. THERMODYNAMICS 



The model is not yet defined in a unique way, because we did not define a Hamiltonian. 
Experimentally, the ground state of the polymer is a compact globule at vanishing force. 
An obvious definition is then H = —hk, with h a positive constant with dimensions of 



energy. With this choice, the probability distribution is a Gibbs distribution, in very good 
approximation [h]]. As a consequence no phase transition can occur. 

The microscopic model with H = —hk does not take the excluded volume effect into 
account. This is a long range effect. It is well known that long range interactions can cause 
phase transitions even in finite one-dimensional systems. So it is appropriate to take the 
following definition for the Hamiltonian 



instead of the more obvious definition mentioned above. With this definition, the ground 
state is a compact globule in absence of an external force. But also the excluded volume 
effect is partly taken into account. This definition is also acceptable form a thermodynamic 
point of view for two reasons. The energy is extensive. It is well known that the thermal 
fluctuations of an extensive quantity (like the end-do-end distance) are proportional to the 
thermal energy in classical systems. 

The Legendre transform of S is the free energy G 



The solution of the set of equations dG/de = 0, dGjd\i = for (3 and F is given in the 
appendix [XI In the thermodynamic limit this solution simplifies to 



The most general case when e ^ /i corresponds with a persistent random walk with drift. 
A persistent random walk is obtained with the choice e = //. This implies F = but non- 
vanishing (3. Also random walk with drift is a special case. This corresponds with e + \i = 1 
and implies (3 = but non- vanishing F. Simple random walk is obtained with e = fi = 1/2. 
In this case both (3 and F equal zero. 

The set of equations (jHJ) or ( 1A2I) has multiple solutions for most values of f3 and F. The 
stable solution is the one which minimizes G. 



H = - 2 (x-(x))\ 
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FIG. 3: (color online). Plot of the force-extension relation at constant temperature, for different 
values of the temperature. The length of the chain is 10 for all three curves. The value of h(5 
equals 0.6, 0.651 and 0.7 for the red, blue and black curve respectively. 

IV. RESULTS 

By inverting the equations (0 [5]) or (El IA2I) one can calculate the force-extension relation 
at constant temperature for any length of the chain. This relation is the central result of 
the present paper. It is shown in figure [U The black solid line corresponds with the stable 
solution (the solution which minimizes the free energy). For small forces, the end-to-end 
distance grows approximately linear with the applied constant force. At Fh/a ~ 0.44, the 
polymer undergoes a transition from a compact globule to a stretched coil. This is clearly 
visible in the figures [1] and [2] because of the sudden increase respectively decrease of the 
average end-to-end distance and average number of kinks. After the transition, the force 
needed to extend the polymer further grows in a strongly non-linear manner. 

The free energy has two minima below the critical temperature. One of these minima 
corresponds with the compact globule, the other with the stretched coil. When the free 
energy of these two states is equal, they coexist. Above the critical temperature, the phase 
transition disappears. This can be seen in figure El where the force-extension relation is 
shown for three different values of the temperature (T < T c , T = T c and T > T c ). 

In figure H] the phase diagram is shown for different lengths of the chain. The critical line 
is first order, up to the critical point. At this point the phase transition is second order. The 
critical line is obtained by numerically solving the set of equations © respectively (IA2I) . 
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FIG. 4: (color online). The phase diagram for different values of the chain length. The critical line 
moves to higher forces when the length of the chain increases. 

Below the critical point, these sets have three solutions. They correspond with two minima 
and one maximum of the free energy. At the critical line we have the extra condition that 
the values of the two minima have to be equal. This is also checked numerically. 

The unstable regime includes two points at which dF/d(x) = (see figure [1]). At the 
critical point these two points coincide and the intermediate regime, where dF/d(x) < 0, 
vanishes. As a consequence, the critical point can be obtained by solving the set of equations 
dF/d(x) = 0, d 2 F/d(x) 2 = 0. This can be done numerically for every length of the chain. 
In the thermodynamic limit one obtains following values for the the critical temperature 
and critical external force T c /h « 1.557 and F c h/a ~ 0.5927 



To summerize, the simple one- dimensional model, introduced here, exhibits a first order 
phase transition. Our approach is a generalization of the mean field approximation. One 
of the methods to perform the mean field approximation starts by assuming a probability 
distribution of the product form, neglecting the correlations between successive units of the 
polymer. The parameters of the probability distribution are then varied to minimize the 
free energy of the total chain. We improve on this scheme by starting with the probability 
distribution of an exactly solvable model. This is an improvement because now the corre- 
lations between successive units are not neglected. The probability distribution depends on 



V. DISCUSSION 



two model parameters e and /i. Next e and \i are varied to minimize the free energy of the 
total chain, as calculated using the correlation functions of the exactly solvable model. 

In this methodology we need two models. At one hand the exactly solvable model whose 
correlation functions are used. At the other hand the more realistic model, which is treated 
in the mean field approach. The exactly solvable model is a random walk with a one step 
memory effect. This model is used in [ll|] to compare the outcome of two different single- 
molecule experiments. In that paper the Hamiltonian is proportional with the number of 
kinks. In the present paper, we increased the physical relevance by defining the Hamiltonian 
by (j3D- In this way the first order phase transition, which is seen in force-clamp experiments, 
is included in the theory. To increase the physical relevance even further we should look for 
mean-field approximations of realistic Hamiltonians. For example, in numerical simula- 
tions are performed with a Hamiltonian which contains a quartic potential in combination 
with a Lennard- Jones potential. The analytic treatment of such a realistic Hamiltonian 
poses technical problems and it is not yet clear how to overcome these. 

The model under study describes a one-dimensional finite system. Therefore, no phase 
transition is expected. However, the present model does exhibit a phase transition. The 
reason is the introduction of a mean-field-like interaction, see the definition ([3]) of the Hamil- 
tonian. It is well known that such interactions can cause phase transitions even in finite, 
nne-d—al Sy8tem , For stance, in flQ a one-di—al, infinite Isin g m „de. is 



considered and the Hamiltonian contains a long-range, mean-field-like coupling. The canon- 
ical phase diagram is derived in and exhibits a phase transition. In [l^] the model 
is further analyzed within the microcanonical ensemble. The authors show that the phase 
diagrams of the canonical and the microcanonical ensemble are not equivalent for this model. 

In force-clamp experiments, one is able to study the mechanical unfolding of polymers in 
the fixed force ensemble. Both experiments 1] and numerical simulations show that a 
certain protein (called ubiquitin) unfolds in a single step. This is in qualitative agreement 
with the results of our one-dimensional model. In order to compare theory and experiment 
quantitatively, our model has to be generalized to higher dimensions. This is a technical 
complication. Our present model cannot be used to discuss the results of experiments 
performed in the fixed stretch ensemble 0, Q]. The reason for this is that the important 
mean-field-like term of the present model becomes zero in this ensemble. Our present model 
can also not produce an unfolding process in multiple steps, which is observed in 



To describe these experiments, the excluded volume effect has to be taken into account 
more accurately. One obtains then a so called self-avoiding walk. In [ijl such self-avoiding 
walks are studied numerically in the fixed force ensemble. The authors show that at low 
temperatures the force-extension curve has a multi-step pattern due to finite size effects. 

APPENDIX A 



In this appendix the results of two straightforward but tedious calculations are given. 
The average of the squared end-to-end distance is obtained with the generating function 
method. The result is 
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The set of equations dG/de = 0, dG/dfi = is solved for j3 and F for arbitrary n. This 
gives 
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